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Abstract 

We investigate the possibility of observing very small amplitude low frequency solar oscillations with the proposed laser 
interferometer space antenna (LISA). For frequencies u below 3 x 10~* Hz the dominant contribution is from the near zone 
time dependent gravitational quadrupole moments associated with the normal modes of oscillation. For frequencies i/ above 
3 X 10~* Hz the dominant contribution is from gravitational radiation generated by the quadrupole oscillations which is larger 
than the Newtonian signal by a factor of the order {2TTrv/c)^ , where r is the distance to the Sun, and c is the velocity of light. 

The low order solar quadrupole pressure and gravity oscillation modes have not yet been detected above the solar background 
by helioseismic velocity and intensity measurements. We show that for frequencies u < 2 x 10^'' Hz, the signal due to solar 
oscillations will have a higher signal to noise ratio in a LISA type space interferometer than in helioseismology measurements. 
Our estimates of the amplitudes needed to give a detectable signal on a LISA type space laser interferometer imply surface 
velocity amplitudes on the sun of the order of 1 — 10 mm/sec in the frequency range 1 x 10~* — 5 x 10"'' Hz. If such modes 
exist with frequencies and amplitudes in this range they could be detected with a LISA type laser interferometer. 
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I. INTRODUCTION 



The proposed ES A/NASA gravitational wave laser interferometric space antenna (LISA) [3,01 will give a unique window 
into gravitational wave physics. LISA will be sensitive to gravitational waves in the frequency range between 10~* and 1 Hz, 
a range currently inaccessible on the ground due to seismic noise. LISA will consist of three spacecraft at the vertices of an 
equilateral triangle of sides 5 x 10^ metres. The system would be maintained in this configuration by arranging that the plane 
of the detectors has an inclination of 60° to the ecliptic and counter rotates with the same period as it orbits the Sun. Any 
two arms will constitute a Michelson type interferometer, a mother spacecraft will send a laser beam to the other two satellites 
where the signal would be coherently transponded back. The interferometer readout would be obtained by interfering the 
incoming signal with the outgoing one and comparing the fractional change in phase shift between the two arms. 

Current studies of the sensitivity of the LISA experiment [3| indicate that the instrumental noise, due dominantly to the 
residual uncompensated accelerations, would be of the order 3x 10~^*^/v'Hz at 10~* Hz (in units of dimensionless strain per root 
Hertz). In addition there is likely to be background "confusion noise" from binary systems which has a comparable magnitude 
at 10"* Hz [3,13, [El- With one year's observation the dimensionless strain that could be detected at a signal to noise ratio of 5 
is estimated to be around 10"'^^ — 10"^". 

The primary goal of LISA is to detect gravitational waves from individual sources (close binary systems, neutron star or black 
hole coalescence) and any stochastic background due to the superposition of waves emitted by binary systems, and possibly 
from the early Universe [1, |3, [3l- Apart from the primary goal, in general, a LISA type interferometer is sensitive to any 
variations in the gravitational field in the frequency range 10~* — 1 Hz. Since the Sun is known to be oscillating in normal 
modes of small amplitude with frequencies in this range, if the amphtudes are large enough the oscillating external gravitational 
field could contribute to the signal detected by LISA [ol. [lol. [Til, [l^ (previously, independently suggested in [H, [3). 

At the present time solar oscillations can only be detected through surface variations in velocity and luminosity. The first 
evidence of surface layer solar oscillations dates back to the work [lEl. The low order gobal oscillations, which are the oscillations 
of interest here, were detected by [l^.[T7| as resolved peaks in the power spectrum of a time series of measurements of the Doppler 
shift of a K and Na line using the integrate light from the Sun. As a result of ground based and space based observational 
programmes upwards of 10^ oscillation modes have been identified in the frequency range 10"'^ — lO"'^ Hz.; knowledge of these 
frequencies has been used to infer the acoustic and dynamical structure of the Sun (pressure, density and rotation as a function 
of radius), placing constraints on the physics of the solar interior and on models of solar (and thereby stellar) evolution. 

The relation between the surface amplitude in velocity and the magnitude of the oscillating gravitational field depends on 
a detailed understanding of the behaviour of the oscillations in the outer layers of the Sun which is not well understood. 
Such modes (g-modes and low order p-modes) are difficult to detect above the solar background noise, which increases at low 
frequencies. Great efi'ort is currently being expended in the search for such modes. In the paper [Q the authors claim to 
have detected a g-mode at a frequency of 1.5 x 10~^ Hz with an amplitude of the order 1 cm, but others have not detected 
any such modes at this amplitude [l9| . Since only quadrupole modes could have amplitudes large enough to be detected by a 
gravitational wave detector at 1 a.u., we confine our study to such modes. 

At a distance of r = 1 a.u. these oscillations are in the near zone (Newtonian regime) for frequencies v such that c/2nv > r 
(i.e. V < Vr ~ 'i X 10~* Hz). The observational implication of the near zone oscillations for LISA was studied in [^. For 
larger frequencies, 1 a.u. is already in the wave zone which means that in addition to the time dependent external Newtonian 
gravitational field the time varying gravitational quadrupole moments also generate gravitational waves, which could give a 
detectable signal on a LISA type interferometer for frequencies v > Vr ~ Z x 10~* Hz [lol. [Til, [l^ . 

In the present work we investigate the possibility of detecting such low frequency quadrupole oscillations with a LISA type 
Iciser interferometer, including both the Newtonian near zone perturbations and the associated gravitational wave emission 
in the frequency range 3 x 10~^ — 10"'^ Hz. We compare the gravitational signals detectable by laser interferometry with 
the velocity signals detectable by whole disc helioseismolgy, and demonstrate that low frequency quadrupole oscillations with 
surface velocity amplitudes below current helioseimic limits could nevertheless be large enough to contribute to the signal 
detected by LISA. If such modes are first detected by helioseismic techniques then the measured frequencies (and predicted 
power) would provide a valuable calibration tool for LISA. On the other hand, if they are not detected by helioseismic means 
then one can look upon LISA as a potential telescope for studying the deep solar interior. 

The plan of this paper is as follows. In Section [Hj we begin with a discussion of the properties of the solar oscillations 
and quantify the relationship between the surface radial amplitude of an oscillation, its quadrupole moment and its horizontal 
amplitude. We express the external gravitational field in terms of the quadrupole moment tensor expressed as a sum over a set 
of basis tensors (corresponding to surface harmonics) which enables us to relate the generation of gravitational waves to the 
surface amplitude of the oscillation. In Section HITl we derive the expected phase shift in an interferometer arm due to both the 
time dependent Newtonian field and to the associated gravitational waves for a given oscillating quadrupole moment. These 
results are then used in Section HVl to determine the response of LISA. Our analysis shows that above frequencies i/~ 3x10"* Hz 
the signal is dominated by gravitational waves. Next, we determine the magnitude of the velocity signal from each mode in 
terms of surface amplitude. In Section|V]in order to analyze the prospects of detectability, we consider the background noise for 
both velocity and gravitational detectors and construct the signal to noise for both gravity and velocity experiments for a given 
assumed frequency resolution. The ratio of these S/N is then independent of the assumed surface amplitude of the oscillation 
and if this ratio is greater than 1 the modes are easier to detect by a gravitational laser interferometer than by helioseismic 
experiments. The outcome of these calculations is summarised in Figure [S] It follows that, for frequencies v <2x 10"'* Hz the 
most of the quadrupole modes are more readily detected in by a LISA type interferometer than by helioseismic experiments. 
Finally, we present our conclusions in Section IVll 
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II. SOLAR OSCILLATIONS AND METRIC PERTURBATION FIELD AROUND THE SUN 



The solar oscillations are normally expressed in terms of a surface harmonic and Fourier time decomposition with any 
variable, for example the radial displacement 5r, expressed in the form 



Sr{r,t) 



E E E a.™(r)S,^(e,0)e'-*, (1) 



n — — oo 



where u; — ujnim are the cyclic eigenfrequencies of modes corresponding to a particular surface harmonic Sim(0, (f>), n labels the 
order of the mode (|n| is essentially the overtone number, the number of nodes in the radial direction), („tm{r) the corresponding 
dimensionless eigenfunctions and (r, 9, (j)) spherical polar coordinates with origin at the centre of the Sun. The modes are 
classified as p (pressure) modes with frequencies increasing with increasing n, and g (gravity) modes with frequencies decreasing 
with increasing n, (for clarity we use negative n for the g- modes in ((TJ). If the basic unperturbed state is spherically symmetric 
the frequencies are independent of azimuthal order m. Rotation lifts this degeneracy giving frequencies ujnem ~ coneo+mQ where 
Cl ~ Qq ~ 3 X 10~® rad/sec is a weighted mean of the solar angular velocity. With a frequency resolution of Av ~ 3 x 10~* Hz, 
as envisaged for LISA (with 1 year's observation time), these individual m- value modes should be resolved provided there is 
sufficient power and the line widths are sufficiently narrow. For frequencies u = a;/27r ~ 10"* Hz, which is the region of interest 
in the present analysis, Qq/lo ~ 10~^ and the eigenfunctions of the modes may be taken to be independent of azimuthal order 
m. The solar rotation axis is inclined at an angle of about 7° to the ecliptic plane, which introduces an additional (small) 
modulation which is neglected in the present analysis. 

There is currently some debate over whether or not any oscillations with frequencies u ~ 10"* Hz have yet been detected 
(c/. [3). The measured line widths at higher frequencies decrease with decreasing frequency and crude extrapolation from the 
measured range suggests Au < 3 x 10^* Hz at frequencies u ~ 10^^ Hz and consequently mode lifetimes more than 1 year. 
We shall assume here that the lines are narrower than the frequency resolution and the modes may effectively be considered as 
monochromatic. Since the external gravitational potential of a multipole of order £ decreases like l/r^"*"^ the modes with £ — 2 
will dominate at 1 a.u.. For this reason we shall only consider the quadrupole modes £ = 2 in the present work. Furthermore, 
since the oscillation velocities are very small compared with the velocity of light only quadrupole gravitational radiation will 
be significant. 

The external Newtonian gravitational potential of the oscillating Sun can then be expressed in the equivalent forms 

f/(r,t) = -G I ^ = f/o(r)-|p"^V.^(i) =t/o(r)-GA/0i?| ^ :^52™(0,0), (2) 

where Uo{r) is the time independent potential, {r,9,<p) are spherical polar coordinates, S2m{d,(j>) surface harmonics of degree 
£ = 2 (normalised to unity over a sphere), Jnm oc e"^^ the dimensionless time dependent quadrupole moments corresponding 
to eigenmodes with cyclical frequencies uj = uj„m- In the above expression V"'^ is the quadrupole moment given by 

/oo 2 
(Sx^x" - S"\^X^'J pit)dV ^ MqRI J2 E CrnJr.mI^^, (3) 
^ n— — oo m — — 2 

where, correspond to a Cartesian coordinate system (with along the 9 = and along (f) = G, 9 = -k 12 axes), the 
Kronecker delta, V^/s = / dx°' dx^ , and are the set of trace-free basis tensors corresponding to surface harmonics and 
Co = —1/5/471, Cm = ^15/47r (for m / 0). We refer the reader to Appendix 1X1 for details. 

The properties of the eigenmodes were computed using the standard Aarhus Solar Model SI [lOl, with surface radial velocity 
amplitudes normalised to 1 cm/sec at the solar surface r = Rq. The amplitude of oscillating quadrupole moments J„. are given 
in Figure [1] (we have suppressed the subscript £ = 2 and m, since the eigensolutions are independent of m). Table |I] gives a 
summary of the results obtained for the amplitudes of relevant quantities derived in this and subsequent sections. Column 1 
gives the radial order n of the mode, column 2 the frequency u, column 3 the horizontal displacement eigenfunction at the solar 
surface C,h, column 4 the quadrupole moment J and column 5 gives the total kinetic energy of the mode E. 

The sun is rotating with a period ~ 27 days so a frame of reference fixed relative to the oscillating sun rotates relative 
to an inertial frame, producing a modulation of the time dependent gravitational field. The solar rotation axis is inclined to 
the ecliptic plane, and hence to the orbit plane of the gravitational detector, producing a further modulation of the signals 
from the oscillations and a very low frequency signal {Av ~ 4 x lO^'' Hz) from the static quadrupole moment induced by the 
rotation. As the inclination is small (about 7°) we neglect these effects in the present analysis taking the rotation axis of the 
Sun perpendicular to the orbit plane of the detector. 

Since the monopole gravitational field of the Sun is weak (GMqItc? ~ 10~* at 1 a.u.), and that of the time dependent 
quadrupole moments and any associated gravitational radiation even weaker, the metric of space-time in the neighbourhood of 
the Earth is adequately described by the weak field limit as 

ds^ = g.^dx' dx^ = (77,fc + h,u)dx^ dx*' = ( 1 + ^ ] <?dt^ ^ ) + h^^dx'^dx'^, (4) 
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FIG. 1: Normalised quadrupole moments J„ (in units of GMRq) for quadrupolar solar oscillations with surface displacement 

UR) = 1- 

where r]ik — diagil, —1, —1, —1) is the Minkowski metric, \hik\ 1, indeces are raised and lowered using the Minkowski tensor, 
and a summation over repeated indices is implied. Roman indeces i, fc = 0, 1, 2, 3 whereas Greek indeces a, /3 = 1, 2, 3. The hik 
have a contribution from the time dependent quadrupole moments 

jv 25{/ GMqRq ^al3^ GMqRq (1\ sr^ sr^ 

h^k = —&ik= SikV ^V^fii-j = SikVc.f){-\ 2^ Cn.Jnn.Tj', (5) 

^ ^■'^ri — — oom — — 2 

and a space-like contribution from the gravitational quadrupole radiation, which with time dependence oc e*^* is 
.Gw 2G <f f~ ^ 2co^G~ 2u^GMqRI ^ ^ 

n— — oo m— — 2 

In the above expression T>ap is the retarded, transverse (i.e. lying in the plane perpendicular to the radial direction of propagation 
of the waves) and trace- free part of Vap. are the retarded, transverse, trace-free tensors corresponding to the basis tensors 

-^m • 

The subdivision in expression Q of the metric field into the Newtonian potential and gravitational wave field at the border 
between the near and the wave zones is obviously a serious over simplification. Nevertheless, this subdivision picks out the 
essential physical aspects of the problem, and there are no obvious reasons for the exact analysis to yield qualitatively different 
answers. In a recent paper [SJ] the authors conducted a thorough analysis of light propagation through the intermediate zone 
between the wave and the near zones. In follo wing papers we hope to analyze the problem of solar oscillations and their 
detectability by LISA using the exact analysis in [2l| |. 

III. PHASE SHIFT IN A PERTURBED GRAVITATIONAL FIELD 

In order to analyze the response of a LISA type interferometer to solar oscillations we firstly need to study the relative phase 
shift for light traveling along the arms of the interferometer in the perturbed gravitational field Q. In order to proceed, let 
us consider an electromagnetic wave of frequency uje and wave 4-vector ki propagating along the arm of a detector from an 
emmitter at coordinate position x\ to a receiver at coordinate position ig. The phase of the wave at is given by 

pB pB pB 

I kidx'' ^ {kiBXg - kiAx\) - x^ dki = kiAixg - xa) + {xg - x^)dki. (7) 
Ja J a J a 

In Minkowski space the wave vector fc* is constant = k\ = LJe(l,n")/c and the phase difference (j>{B) — <j>{A) is simply 
kiA(x^B ~ ^\)- The second term on the right hand side Q is identically zero in the Minkowski limit. 

In the time-dependent space time Qik the first order perturbation in phase, due both to the departure from Minkowski space 
time and the displacement of the receiver B (taking A as fixed) is given by 



S<4i ^ kiAS{x''B — x^) + I {x'b — x')dki 

J A 



(8) 
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TABLE I: The frequency u (Hz) and the respective amplitudes of horizontal surface displacement Ch{Ro), quadrupole moment 
J, the total kinetic energy of the mode E (gm ■ cm'^/sec^), gravitational signals Sm and velocity amplitudes Vm (m/sec), for 
unit radial surface displacement (^^(-Rg) = 1 as a function of the radial order of the mode n (negative and positive n correspond 
to g-modes and p-modes respectively). The last three columns show the expected instrumental noise Bi (l/\/Hz), binary 
confusion noise Bb (l/\/Hz) and surface velocity noise B„ (m/sec/VHz) respectively for the LISA mission at the corresponding 
frequencies. 



We take A as the origin of coordinates so x\ = 0, and take the unperturbed ray to be given by x' = Afc*, where A is an afHne 
parameter varying from at A, to A at B. Since = dx^/dX, the null geodesic equation in this weak field approximation 
reduces to 



dk' 



and the last term in can be expressed as 

/ {A-\)k'dh = ^ (A - A)/i™,ifc'"fc"fcVA = ^ / (A- A)^^fc'"fc"dA. 
Ja 2 Jo 2 Jo dX 

Substituting this result into ([8]), defining P — (xg — x^j^, 5P = S{x^b — xa), and integrating by parts gives 



Sip — kiASt -\- 

SF is given by the geodesic deviation equation 



-(A - X)hmnk'^k'" 



A ^ , 



k"'k"dX. 



(9) 
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which, in the slow motion approximation appropriate to the current analysis, reduces to 



^^((5? ) = Rokol 



Since in the weak field approximation 



(10) 



(11) 



From jni and it follows that R°^^ = and (fSr/dt'^ = 0, respectively. Recalling that hoa = 0, and P = Tq + SP with 
lo = (l,n'*) we obtain 



df2 



r 



where V^^ = d'^ / dx^ dx'^' . In view of the time dependence of hik oc e'"', this equation can be integrated to give 



Now since k" — ujeu" /c, A = cl/uje, where is the frequency of the electromagnetic wave, substitution into ([9]) gives 



2c V ^ ' J 2c 2 Jo 



(12) 



(13) 



(14) 



where we have introduced I = Ualo as the unperturbed distance between A and B. 

We now write S(fi = S^n + Sificw, where 5ipN is the contributions from the Newtonian potential 5U and Sifcw is the 
contribution from gravitational waves. Using expression (|4]) and expanding the hmn in (|14|) by Taylor series yields 



(15) 



Taking into account ((5]) and ([6]), retaining just the leading terms in powers of uol/c <^ 1 and l/r <^ 1 in (|15[) . the phase shift is 
given by 



= n n 

2c 



2 

-5- Vaflftoo — ft. 



\ _ ijj G 
c 6lj^ 



(16) 



where = / dx" dx^ dx'' dx^ . Since the unperturbed phase shift is ^ = uJe.l/c, the fractional change in phase shift can be 

expressed in the form 



G 



Sip 

ip J 2ijp-r'^ 



a0 



(17) 



where Nafi = rianp, and 

Tap = i N^^^afj-ft = 2nanp — 10ft(7ian,3 + ncTi/j) + 5(7^^ — l)nan,3, (18) 

and n" is the unit vector in the radial direction and ^ = n^n^. The details of the derivation of Tap are given in Appendix iBl 



IV. RESPONSE OF THE INTERFEROMETER 



A laser interferometer detector such as LISA consists of 3 arms AB, AC, CB, in circular orbit around the Sun. The 
interferometer's response is given by the difference in the fractional change in phase shifts of the round trip signals along any 
two of the arms, for example ABA and ACA, as 

where A indicates the difference between the arms AB and AC. 

The response 5 is a function of the angles Xi'^i'/'j^i Ihe quadrupole moments Jm and the frequencu lj. x is the angle 
between the radius vector to the detector and the instantaneous orientation of the reference axes corotating with the sun, l 
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the inclinations of the plane of the detector to the orbit plane, ^ the angle between the arm AB and the direction of the orbit 
and V the angle between the two arms of the detector. In the LISA concept, l = v = 7r/3 and the detector rotates with the 
same period as it orbits the Sun. The angle tp decreases as the position angle of the detector (on its orbit relative to a fixed 
inertial frame) increases, i.e. ^ = tpo — (j) where tpo is a constant equal to the orientation of the detector relative to the orbit 
plane at the arbitrary zero of the orbit angle <p. 

We now express the quadrupole moment tensor in terms of the trace- free basis tensors J'^^^ (see (|3])) and determine the 
response 5*^71 for both gravitational waves and the Newtonian signal for quadrupole modes of azimuthal order m and frequency 
v — Lu/2n. Introducing the effective gravitational amplitudes 

we have 

where the Newtonian and gravitational wave contributions to the signal can be represented in the form 



'e// -3^mJm (■^) (-^) (l^) /m^i 

(20) 



where (see Appendix [CI for details) 
and 



"-i^ff - ^mJm \ —r) \ — j Im, 



rN 7"CK/3 A 7^ 



(21) 



10~''Hz, = ( —) 3 X lO'^Hz. 



27r\/ ' \2-KrJ 



Note that, as follows from (|20p . for and of comparable order the gravitational wave contribution dominates for 

> i/,. ~ 3 X 10"'' Hz. 

To determine the functions , fm^' consider the detector to be at the point P on its orbit with orbit angle (j!> relative to a 
fixed inertial frame, and let x be the angle between the radius vector to P and the axis of the corotating system x°'. We 
then define local transverse Cartesian coordinates at P such that is in the outward radial direction, in the direction of 
the orbit and ^"^ perpendicular to the orbit plane. The basis tensors in the coordinates are given by 

where 

= (cosx,sinx,0), = (- sin x, cos x, 0), = (0,0,1), 

are the unit vectors along the axes in the x" coordinate system. The transverse trace-free tensors are then given by (see 
Appendix [C} 

f23 _ 732 _ r23 r22 _ _ ?33 _ }_ ( t22 _ t33\ fl-y _ 
-'m -'m 5 m 2 V /' ^ 

In the coordinates, in light of the definitions of angles t, tp and v, the unit vectors ng,nj are given by 
= (— sin tp cos i, cos ip, sin-;/; sin t), = (— sin ip' cos i, cos ip' ,smip' sin t), ip' = ip + v. 
In light of the definition the coordinate system the unit radial vector ff is given by 

n-' = (1,0,0), 

After some laborious algebra we determine the tensors AN-y,,, AT-y,, in the coordinates, and hence the functions f^ifrn^ 
defined in (|21|l . The solutions for arbitrary L,v,ip,x are given in Appendix [Cl 
For LISA, setting l — v = n/3 and tpo — (i.e. tp = —(p), we obtain 

/iv = 3^ 3in(2^ + 7V /3), /o«^" = sin(2V + ^/3), 

/f = -3 [2 cosx sin(2V'-f 7r/3) -sinx cos(2i/> -f 7r/3) ] , fi^ = § sinx cos(2V' -f 7r/3), (22) 
fi^ = [f cos2x sm{2iP + 7r/3) - 8 sin2x cos(2V' + vr/S)] , f2^' = cos2x sin(2i/; + vr/S). 
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The values for m = —1, —2 are obtained from those for m = 1, 2 by replacing X by X ^ it /2m. 

Relative to a fixed inertial frame at time t the detector is at orbit angle = and the corotating axis is at an angle 
^Qt, where Q.d is the angular velocity of the detector around the sun and the angular velocity of the Sun relative to this 
inertial frame. The angle between the radius vector from the sun to the detector and the rotating axis is therefore x = ^at 
where Q.s = 2-n /Ps where Ps ~ 26.75 days is the synodic period of solar rotation, that is the period relative to a reference frame 
orbiting the sun at 1 a.u. Thus, substituting t/) = —<j) and x = {^a/^d)4' into (|22[) gives fm{4>) fm^ ((f)) as a function 
of the position of the detector on its orbit determined by the orbit phase angle 0. In Figure [J] for illustration, we show the 
fm{4>) smd the fm^{4>) modes. Since m = correspond to axially symmetric modes, the functions /o(<^) are independent of 
the rotation of the Sun. The other modes (m — ±1, ±2) display the modulation of the signal due to solar rotation. Combining 
the contributions from and f'^^ gives 



1 I 7 



i(27/; + ^/3). 



[ + I rr) cos(2V' + tt/3) - 2 cos x sin(2V' + 7r/3)] 



(23) 



where 



S2iu,^) = 2V3C2 J2 # ^ [sin2x cos(2V' + tv/3) - § (l + X cos2x sin(2V' + ^/3)] , 



-(t>, x = <t>Pd/PG~lS.65cl,, Co 



Cm — 



m = ±1,±2. 



Again, in a similar fashion to (|22p . values for m = —1, —2 are obtained from those for m = 1, 2 by replacing x by x — 7r/2m. The 
expressions in (|23|l contain both the contributions, from the Newtonian potential and gravitational waves (terms proportional 
to {u/vr)^). As anticipated above, the gravitational wave contribution dominates when > i/^ ~ 3 x 10~* Hz. 

The amplitude of the signal of given (m, v) varies with the rotation of the sun around its axis and over the coarse of the 
year. For an observation time T ~ 1 year the amplitude is given by the root mean square averaged over T. The resulting 
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FIG. 3: Limits on detectable signal from LISA with signal to noise S/N=5. (a) Instrumental noise, (b) binary confusion noise 
(Bender and Hils 1997). 



amplitudes are given in Table |T] (columns 5, 6 and 7 for m = 0, 1 and 2 respectively) for values of the quadrupole moments 
given in column 4. These signals scale linearly with the value of the quadrupole moments. 



V. DETECTABILITY OF GRAVITATIONAL SIGNALS FROM THE SUN AND COMPARISON WITH 
VELOCITY EXPERIMENTS 

The expected sensitivity of the LISA experiment for a single 2-arm detector has been evaluated in [22| . At low frequencies 
the dominant contribution is from the acceleration noise cx at intermediate frequencies from shot noise, and at high 

frequencies the sensitivity declines when the path length along the detector becomes comparable to and greater than the 
wavelength of the gravitational wave. In the frequency range of interest this instrumental noise level is given by 



/i/(!^)~^ 2x10"*' +1x10""'— NHz, (24) 
\smi/ \ I 

where i = 7r/3 is the angle between the arms of the interferometer. In addition to the instrumental noise there is also a 
"confusion noise", hsiy), due to the integrated gravitatonal wave contribution from galactic and extragalactic binary systems. 
The magnitude of this "confusion noise", /is(!/), is uncertain but has been estimated in [3| and their estimate is given in column 
12 of Table 1 and shown in Figure O 

Thus, for LISA, the threshold of detectibility hTD at a signal to noise level SjN , when direction to the source is known, is 
given by 

hTD^{h,+hB)(^{^. (25) 



^fT) \N 

where T is the time of observation (in sees). 

Let us now turn to the sensitivity of velocity experiments. The solar oscillations of low degree are obtained by taking the 
power spectrum of a time series of measurements of the Doppler shift (or velocity) of a given spectral line in the integrated light 
from the sun. The frequencies are the peaks in this power spectrum. The velocity amplitude of a given quadrupole oscillation 
is therefore given by integrating the component of the surface velocity in the direction of the observer over the visible solar 
disc, 

where k" is the unit vector in the direction of the observer, ^ — k"fa is the cosine of the angle between the unit radius 
vector fa and the direction to the observer, and /(^i) the appropriate solar limb darkening function which incorporates the 
angular dependence of the intensity of radiation at the solar surface (see for example [l^). We here take /(/i) = a + with 
a — 0.55, b — 0.45, which is a reasonable approximation for the whole disc velocity measurements by the GOLF [Sj and BiSON 
[25j experiments. 
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FIG. 4: Solar background noise in velocity (m^/sec^/Hz) as determined from the GOLF experiment on SOHO. 



In a system of spherical coordinates (r, 9, cj)') corotating with the sun, the velocity on the solar surface due to an oscillation 
mode is 

where C,r,C,h are the dimensionless radial and horizontal displacement eigenfunction at the solar surface R — Rq. Let the 
corotating axes be at an angle 00 to the direction to the observer and take cj> = cj>' — 4>i.) then 

k°'va = Vr sin 9 cos (j) -'rVg cos 9 cos 4> — "Vtf, sin (j>, fi = sin 9 cos (j)- 

Evaluating the integral (|26|) gives the velocity amplitudes 



/5 /^ (16a+15b)Cr + 6(8a + 5b)Ch ^ /- / cos2,^o A 

For m — ±l,S2m oc cos 9 is antisymmetric about 9 — 71/2 and the integrals p6p for ^±1 are identically zero. The m — ±1 
modes are not detectable in integrated velocity in the equatorial plane. 

The surface horizontal displacement is known in terms of the surface value of the radial displacement eigenfunction from 
integration of the eigenvalue equations for the oscillation (column 3 in Table |T]). Since the normalised value of the quadrupole 
moment is J2 <C 1 for all modes (column 4), it follows from the equations governing the oscillations that ~ Cr g/^'^ Rq (see 
[2^). The velocity amplitudes with ~ 1 are listed in columns 9 and 10. The background noise in velocity experiments is 
predominantly from velocity fields on the solar surface (active regions, granulation, meso-granulation, supergranulation) and 
the cumulative effect of these motions has been estimated in [SJ. However in contrast to the gravitational case this background 
noise has been measured in a number of helioseismology experiments and in Figure |4] we show the background velocity noise 
determined by the GOLF experiment on SOHO [13, [2^. This is in reasonable agreement with the value obtained by the ground 
based whole disc networks BiSON [2^, that obtained 900 cm/sec/^/Hz at u — 5 x 10^** Hz and IRIS experiment (s^] which 
obtained 2 x 10'^ cm/sec/\/Hz at = 10"** Hz. These values are a factor ~ 5 below the model predictions of In the region 
of interest this background noise Bv {v) can be approximated by 



We now make the following comparison. We assume a frequency resolution in the power series analysis of 3 x 10^** Hz 
(corresponding to 1 year's observation time) and the value of a velocity signal from an oscillation mode that has S/N — 3 using 
the above background noise estimate (|28|l . This gives a value of the dimensionless radial eigenfunction through (|27[) for the 
velocity amplitude for each frequency and m = 0, 2. This then determines the amplitude of the quadrupole moment J and the 
gravitational signal strength S obtained by simply scaling Sm in Table |T] by the value of C^r- We then compare this gravitational 
signal strength with the background noise for the LISA experiment computed from (|24|l under the same assumptions. The 
results are shown in Figure [5] We see that for frequencies < 2 x 10"'* Hz the signal from the m=2 modes are stronger in the 
gravitational experiment than in the velocity experiment. In Figure [5] we also show the S/N for the gravitational experiment 
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FIG. 5: Comparison of signal to noise (S/N) in LISA experiment for S/N = 3 in solar velocity experiments, for oscillation 
modes with = 2, m = 2. Solid line shows the S/N when only the instrumental noise is included, while the dashed line shows 
the S/N when both the instrumental noise and the noise from binary systems is included. 



when binary confusion noise is included. Note that since the averaged velocity for the m = ±1 modes is zero there is no limit 
on the gravitational signal set by the helioseismology experiments. 

This comparison between the 2 experiments is of course independent of the assumed S/N in the velocity experiment, what 
is being compared is the detectability of an oscillation mode by the two techniques and is essentially the ratio of the S/N for 
the gravitational experiment to S/N for the velocity experiment. 

We can reverse the comparison and ask what would be the velocity amplitude of modes that are at the margin of detectibility 
with the gravitational detector, and how does this compare with the solar background noise? This is done in Figure |6l here we 
see that the velocity amplitude at low frequencies for the m = 2 mode is well below the noise level of velocity experiments for 
a frequency resolution of 30 nHz. For high frequencies > 3 x 10^'* Hz the gravitational signal is predominantly determined 
by the gravitational wave component. For comparison, the figure also shows values for the velocity amplitude obtained from 
These values were calculated using expression (3) and the data in Table 1 in [^, note that these values are averaged 
over m-modes. As can be seen, the estimates from [y| are in a reasonable agreement (taking into account the averaging over 
m-modes) with our results at lower frequencies ^ 3 x 10~* Hz. Let us compare the results for frequencies <; 3 x 10^^ Hz. 
As was mentioned above, at these frequencies the major contribution to the signal comes from the gravitational waves. This 
gravitational wave contribution was not analyzed in [g]. For this reason, as can be seen on Figure [6] at higher frequencies 0] 
give higher values for the detectability threshold in comparison with present results. In other words, the incorporation of the 
gravitational wave contribution makes detection of solar oscillations by LISA more feasible. 

In the above considerations we have assumed that the signal is monochromatic, or rather that the line is narrower than the 
frequency resolution. This is the prediction from extrapolating the p-modes observed by helioseismology at higher frequencies. 
Since the rotational splitting of the modes is known, signal enhancement techniques can be used, superposing power in frequency 
bins separated by mfi©, as is being done in the search for g-modes in helioseismology by the Phoebus group (sj]. If the modes 
have a line width in excess of the frequency resolution the detectibility in both velocity and gravity is correspondingly reduced 
but the ratio of S/N for the two experiments remains the same. 



VI. DISCUSSION AND CONCLUSIONS 

In this work we have analyzed the prospects of detecting the gravitational signal from solar oscillations in a LISA type 
interferometer, and compared them with capabilities of velocity experiments. At low frequencies, < 3 x 10^'' Hz, the 
quadrupole oscillations of the Sun might be detected in LISA through their contribution to the time varying Newtonian (near 
zone) gravitational potential. At higher frequencies, ly > 3 x 10~* Hz, LISA might observe a gravitational wave (far zone) 
signal from the solar oscillations. For frequencies v < 2 x 10^'' Hz the signal will have a higher S/N in a LISA type space 
interferometer than in helioseismology measurements, while for higher frequencies the signal would be more readily observable 
through helioseismology measurements. 

Although, as was mentioned at the end of Section [Hi our analysis was based on the subdivision of the gravitational field 
into Newtonian potential and gravitational wave field at the border between the near and the wave zones, there is no reason to 
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FIG. 6: Minimum surface velocity amplitude for signal to be detectable by LISA with S/N = 1 for modes with m — (solid 
line) and m = 2 (dashed line) . The thin dotted line show the velocity amplitudes evaluated from . 



believe that an exact analysis will give qualitatively different results. In particular, the characteristic features in signal to noise 
ratio curve on Figure [5] should remain, although there might be a slight change in the hight of the curve around frequencies 
1/ 3 X 10"* Hz corresponding to the intermediate zone. In a following paper we hope to analyze this question in more detail. 

The low frequency solar oscillation modes, considered in this paper, are more sensitive to the structure of the solar core 
than the higher frequency p-modes. Helioseismology experiments have so far been able to measure only the higher frequency 
p-modes. Hence, if the low frequency modes are detected by LISA, they would advance our understanding of the structure 
of this central core. This is because, the helioseismology experiments measure the surface velocity which is dependent on the 
detailed structure of the outer layers of the sun, and this has to be subtraced off to give diagnostic information on the solar 
interior. On the other hand, the interferometer experiments measure the actual quadrupole moments which are determined 
by the oscillations in the deep high density solar interior and therefore encode information about the solar interior which is 
independent of the structure and physics of the outer layers. 

The predicted amplitudes of these low n modes are still uncertain. The p-modes are thought to be stochastically excited by 
the convective motions, and if this is the excitation mechanism of the g-modes the predicted amplitudes of the low frequency 
modes are very small. But as was first pointed out in [13, H^, the steep gradient of ^He in the inner half of the sun could also 
provide an excitation mechanism for these modes which may be damped by parametric resonance with other modes or by mild 
turbulent diffusion, but still be of sufficient amplitude to give rise to a detectable gravitational signal. 

We note also that the frequency of 10"'* Hz corresponds to a period of the order of 160 minutes. There have been repeated 
suggestions that such a signal has been seen in helioseismology experiments, many but not all of the claims being withdrawn 
following more detailed analysis (see [l^l). An analysis of the GOLF data [s^l placed an upper limit on the velocity amplitude 
of such a mode at about 1 cm/sec, an estimate compatible with searches by the Phoebus group [sJl. We note that for such 
a 160 minute mode an S/N = 1 in the helioseismic velocity experiments would correspond to an S/N = 20 in a LISA type 
experiment. 

The search for g-modes in the helioseismology experiments is ongoing and hopefully will result in the confirmed detection 
of some modes prior to the launch of LISA. Were such modes to have been identified, they would provide a valuable known 
signal for the space interferometer which could then be used to calibrate and test the experiment, giving more credence to the 
interpretation of signals from more distant astrophysical sources. 
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APPENDIX A: THE QUADRUPOLE BASIS TENSORS X'^^ AND SURFACE HARMONICS 52m 

The quadrupole gravitational potential can be expressed in the two equivalent forms (see ([5J) 

gMqrIY^ ^-^s^^(e,4>) = G\v-^v^p(S\, (Al) 

m— — 2 n— — oo 

where {r,6,(j)) are spherical polar coordinates, x"" = {x,y,z) Cartesian coordinates with = the z axis and (f) = Q, 9 = 1^/2 
the X axis, and Va/3 = d'^ / dx" dx^ . Jnm are the (dimensionless) quadrupole moments, T)"^ the quadrupole moment tensor and 
S2m are the real surface harmonics normalised to unity over the sphere 

3 cos^ 61-1 / 15 . „ „ / cos \ „ / 15 . 2 „ / cos 20 

82,0 = \ - , S2,±i = a/ — ^ sm26l ■ I , S2,±2 = W 77^ sm e . I, 

V 47r 2 V 167r \^ sm0 y ' V 167r \ sm2(/) 

We define a set of 5 independent trace free tensors as 

/100\ /00l\ /000\ /100\ /OlO 

To = 1 , Ji = , T_i = 1 , T2 = -1 , J_2 = 1 
\00-2/ \100/ \010/ \000/ \000 

Taking into account 



l\ ( ^\ (iXaXp — r^S, 



aj3 



where 5^/3 is the Kronecker delta, and 

x^ = X = r sin 6 cos 0, x^ = y — r sin 9 sin 4>, x^ — z — r cos 9, 

the TJ^'' satisfy the relations 

lo" v.;, (i) = ^{x^ + y^- 2z^) = ^ i^if^ = -I, A/f52,o, 



TfVo.p{^) = ^2xz = ^ sinecos0cos0= ^^^52,1, 

I-? (i) = ^ 2yz = ^ sin^cos^sin,^ = ^ a/#52,-i, (A2) 



^cp (i) = ^ (2;' - t/') = 4 Sin^ eCOs20 = 3 S2,2, 



V.^ (i) = ^ 22^2/ = sin^ esin20 = S2.-2. 



We now expand the quadrupole moment tensor in the form 



2?" — MqRq y ^ y ^ Cm Jnm^Z 



n — — 00 m— — 2 



Inserting the above expansion into the right side of ()A1|) . taking into account (|A2|) , , we arrive at the following expression for 
the coefficients Cm 

^° = -\/5' ^- = \/l --±1'±2. (A3) 

APPENDIX B: DERIVATION OF THE DETECTOR TENSOR T,^ 

We use a cartesian coordinate system x" with x^Xa ~ ■ Then 

V^. (i) =-^+3^, 

V^,.(i)=+3^ + 3^+3^-15ii^, (Bl) 

^i/37e (f ) = ;3-(^E7^/3q + ^tfiS-ia + <5sq(5t,^) 
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Defining the unit vector in the radial direction 

V r r r / 



we get 



.4 / 1\ 1 



Defining ji = rffit and contracting witli rfrC gives 

— 5 (fi/gnQ + ^npfia + ^.nanfB + fmpna + iirianp + ^^Saii) + 35 ji^Ufsn^ . (B2) 

Since this tensor is to be contracted with the symmetric trace free tensor T>^°' = T>°'^ and (5/3qC" = we only require the 
trace free components of tliis tensor. For this reason the terms with 5^/3 can be removed in (|B2|l leaving 



r 



2n0na — lOfiriisria — lOfiriaTif^ + 5(7^^ — l)nf3na. (B3) 



trace— free 



APPENDIX C: THE FUNCTIONS /f (<;/>) , /f ■^(0) 

Here we compute the source tensor in local transverse Cartesian coordinates ^" at a point P on the orbit of the detector 
where (,^,£,^ are in the orbit plane, in the the outward radial direction, is in the direction of motion and perpendicular 
to the orbit plane, x is the angle between and the x^ direction of the coordinate system x" in which the multipole moments 
are determined. The unit vectors along the axes in the x" coordinate system are 

65 = (cosx,sinx,0), e| = (- sin x, cos x, 0), e| = (0,0,1). 

Hence the basis quadrupole tensors in the ^ coordinate system are given by = e^e^T"'' and are explicitly given by 

cosx \ f ^ sinx \ 

-sinx , /-I = cosx , 
/ y sin X cos X / 

(cos 2x — sin 2x \ / sin 2x cos 2x \ 

— sin2x — cos2x , 1-2= cos2x — sin2x . 
00/ \ 00/ 

The transverse-trace-free radiation tensors in the ^ coordinates are then 



















1 


• 






Vo 











000\ /000\ /ooo 

1 , /i = -sinx 1 , /-I = cosx 1 
00-1/ \oio/ VOlO 



/000\ /ooo\ 

72 = -icos2x 1 , J_2 = -isin2x 1 . 

\00-l/ \00-l/ 

Now let the plane of the detector be at an inclination angle t to the plane of the orbit (t = 7r/3 in the LISA experiment), and 
let the detector arm AB be at an angle t/j to the direction of the orbit (^^). In the the ^ coordinate system the the unit vector 
along the arm of the arm AB of the detector is 



= (— sin cos t, cos sin ?/> sin t), 



so the projection tensor A*"^^ = n°n'' is 



Csin^ -0 cos^ L — sin -0 cos cos t — sin^ V sin t cos t 
— sin 7/; cos ?/> cos t cos^ V cos '0 sin -0 sin t 

— sin^ il) sin t cos t cos -0 sin sin i sin^ sin^ i 
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Let the angle between the arms AC and AB of the detector be v so that AC makes an angle + v with the ^ direction, the 



projection tensor N^c is given by the above result with -0 replaced by ■0 + so the difference 



(— sin(2i/; + u) cos^ t cos(2i/j + u) cos t sin(2'i/' + sin t cos i \ 
cos(2?/j + u) cos t sm{2ip + v) — cos{2ip + v) sin l 

sin(27/; + d) sin t cos t — cos(2'i/' + sin t — sin(2?/; + u) sin^ t / 

The signal functions for gravitational waves fj'^ are given by 

f^^' = AN^pi^^ = sin V (am{2ij + v) Tf^ - 2 cos(2V' + v) sin l - sin(2V' + v) sin^ l I^'^ , 

which gives 

fo^ = 1.5 sinu(l + sin^t) sin(2V' + u), 
jGW _ 2 giji ^ sin t sin cos(2'(/) + v), 

jGW _ _2 sijju sint cosx cos(2i/) + u), (CI) 
jGW ^ _Q 5 smv{l + sin^ t) cos2x sin(27/; + v), 
jGW ^ _Q 5 sin„ ^ sij^a gjj^2x sin(2V' + v)). 
The Newtonian signal functions = AT^pIm^ where 

ATc/s = 2AA/-^/3 + WAMc,0 + 5ARci3, 

with 

A^Q^ = iJ.{n''nP + n^n°'), Rap ~ — l)n°'n^ , ^ — fCn^. 
In the 5 coordinates n° = (1,0,0), and for the detector arm AB jj, — n"na ~ — sin i/; cost so 

(2 sin^ ip cos^ t — sin ip cos -0 cos t — sin^ -0 sin t cos t \ 
— sin ip cos t/j cos i I , 

— sin^ ?/> sin t cos t / 



and hence 



(— 2 sin(2-(/) + v) cos^ I cos(27/) + w) cos t sin(27/) + u) sin t cos t 
cos(2-i/) + u) cos t 

sin(2?/' + 1)) sin t cos t 



The tensor Rab = (7/x^ - l)n"n'' gives 



/ 1 

Ai?= -7sinwsin(2V' + u)cos^i 

V 



So finally the tensor ATap is 



C— 17sin(2?/' + w) cos^ t — 8 cos(2?/' + u) cos t — 8 sin(2i/; + t;) sin t cos t \ 
~8 cos(20 + u) cos t 2sin{2ip + v) — 2 cos(27/; + u) sin i I . 

— 8 sin(2i/' + w) sin t cos t — 2 cos(2i/' + u) sin t — 2 sin(27/' + u) sin^ t J 

The signal functions = I^^AT^p are then 

fo' = sin u (6 — 21 cos^ i) sin(2'!/' + ti), 

J]*^ = —4 sin w sin i [4 cos t cos x sin(2?/> + u) — sin x cos{2ip + v)] , 

/il'i = —4 sin u sin t [4 cos t sin x sin(2i/) + u) + cos x cos(2'0 + v)] , (C2) 

= — sin V [(17cos^ t + 2) cos 2x sin(2?/' + u) — 16 cos t sin 2x cos(2?/' + w)] , 
/i^2 = — sin V [(17cos^ t + 2) sin 2x sin(2i/' + u) + 16 cos t cos 2x cos(2i/' + u)] . 
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Setting L — V = n/3 in (|Clfl and (|C2|) brings us to expression (|22|l for the specific case of LISA. 
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